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ABSTRACT 

This paper presents details of a flutter and stability 
analysis of aerospace structures such as hypersonic 
vehicles. Both structural and aerodynamic domains are 
discretized by the common finite element technique. A 
vibration analysis is first performed by the STARS code 
employing a block Lanczos solution scheme. This is 
followed by the generation of a linear aerodynamic grid for 
subsequent linear flutter analysis within subsonic and 
supersonic regimes of the flight envelope; the doublet 
lattice and constant pressure techniques are employed to 
generate the unsteady aerodynamic forces. Flutter analysis 
is then performed for several representative flight points. 
The nonlinear flutter solution is effected by first 
impl e men t ing a CFD solution of the entire vehicle. Thus, 
a 3-D unstructured grid for the entire flow domain is 
generated by a moving front technique. A finite element 
Euler solution is then implemented employing a quasi- 
implicit as well as an explicit solution scheme. A novel 
multidisciplinary analysis is next effected that employs 
modal and aerodynamic data to yield aerodynamic damping 
characteristics. Such analyses are performed for a number 
of flight points to yield a large set of pertinent data that 
define flight flutter characteristics of the vehicle. This 
paper outlines the finite-element-based integrated analysis 
procedures in detail, which is followed by the results of 
numerical analyses of flight flutter simulation. 

INTRODUCTION 

The accurate prediction of flight characteristics of 
hypersonic vehicles is vital during its design stage as well 
as prior to flight testing to ensure flight safety. Such a 
vehicle is expected to exhibit unprecedented levels of 
interaction among various disciplines such as structures, 
aerodynamics, and controls engineering, among others. 
For complex configurations, it is necessary to implement 
an efficient discretization procedure for effective and 
accurate idealization of the vehicle. In this connection, the 
finite element method proves to be a viable technique to 
model both solids and fluids continua and therefore is a 
natural choice for aeroelastic analysis that involves 
interaction of associated disciplines. 


Multidisciplinary research at NASA Dryden Flight 
Research Facility is aimed at developing integrated 
aero elastic and aeroservoelastic (ASE) analysis capabilities 
that can be conveniently employed for effective prediction 
of flight-critical, dynamic stability and control performance 
parameters and characteristics. A unified, nonlinear, 
multidisciplinary simulation analysis approach has been 
developed to accommodate these requirements, thus 
yielding a comprehensive analysis framework. The 
primary elements of this approach include 3-D unstructured 
fluids grid generation with adaptive mesh capabilities, 
finite element structural and heat transfer analysis 
procedures with advanced materials, as well as integrated 
aeroelastic and ASE analysis capabilities. Since such 
nonlinear analyses require extensive computing resources, 
much emphasis has been placed on developing novel 
solution schemes that reduce such solution effort 
significantly. Further, dedicated parallel processors have 
also been used to achieve the desired solution for practical 
problems within a reasonable span of time. 

NASA’s finite-element-based, multidisciplinary 
modeling and simulation computer program, STARS, 1 
developed for the solution of practical problems associated 
with NASA-critical flight projects, is described here in 
some detail. Extensive graphics-oriented pre- and 
postprocessing capabilities render the code useful for the 
solution of complex problems. The associated numerical 
formulations pertaining to various individual disciplines as 
well as the integrated solution techniques are discussed in 
the paper. 

Some numerical examples are also presented to 
demonstrate the accuracy and efficacy of the NASA STARS 
code. The first example relates to a rectangular panel with 
clamped edges, and the STARS nonlinear flutter solution is 
compared with that obtained by linear analysis as well as 
wind-tunnel tests. A complete generic hypersonic vehicle 
is chosen next as an example problem. Relevant details of 
the structural modeling as well as a complete vibration 
analysis are presented, which constitute the first important 
step toward achieving an aeroelastic stability solution. 

The results are presented for a set of parallel aeroelastic 
analyses pertaining to low Mach numbers within the flight 
flutter envelope using linear, unsteady aerodynamic theory 
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based on panel methods. This is followed by a number of 
CFD analyses pertaining to some representative flight 
conditions within the specified flight flutter envelope. 
These results are used in an attempt to generate unsteady 
aerodynamic forces and subsequently the damping 
characteristics that are indicative of aeroelastic instability. 


K,K*,Kg = elastic, centrifugal, and geometrical 
stiffness matrices 

u = unknown displacement vector 

g = structural damping 


DESCRIPTION OF NUMERIC AL TECHNIQUES 


imaginary number. 


VT 


To s inuiiam nonlinear performance characteristics of 
advanced aerospace vehicles, it is necessary to integrate a 
number of relevant disciplines in a consistent fashion. A 
common employment of the finite element technique for 
the various disciplines such as the fluids and solids 
continua ensures accurate simulation of their interactions. 
The various finite element modules of the STARS 
program, developed in this connection, are shown in figure 
1, and some relevant details of their formulations are 
presented next 



Fig. 1 STARS analysis modules. 


f(t) = externally-applied, time-dependent, forcing 
function 

The various matrices in equation (1) are usually large but 
sparse in form for complex practical problems, and the 
STARS numerical analysis module is designed to solve 
such matrix simultaneous equations and eigenvalue 
problems in an efficient fashion. 

Layered, anisotropic, advanced composite finite 
elements are used for the design of advanced aerospacecraft. 
Related data for temperature-dependent properties are stored 
in the material module of the program. 

The dynamic behavior of a viscous, heat-conducting, 
compressible fluid obeying conservation of mass, 
momentum, and energy may be expressed by a set of 
partial differential equations 

3t + dxi = fb * 1,2,3 (2) 

where the solution, flux, and body forces column vectors as 
well as the viscous stress tensor are defined as 

V = { p puj pE } (3) 


Fi = { pui pujUj+p8ij+<Jij 

gj 

Ui(p + pE)+ui<Jji+ k^ 7 } (4) 


Thus, the structural analysis module is capable of 
performing static, vibration, buckling, and dynamic 
response analyses of complex, practical structures having 
general, anisotropic material properties. A typical matrix 
formulation encompassing a broad class of structural 
problems may be written as 

Mu + (C c + C d )u + [K(l+/g) + K' + K g ]u = f(l) (1) 
where 

M = inertia matrix 

C c .Cd = Coriolis and viscous damping matrices 


fb = { 0 fbj uifbi ) (5) 



in which p, p, and E are the density, average pressure 
intensity, and total energy, respectively; Sjj is the 
Kronecker delta; uj is the velocity component in the 
direction xj of a Cartesian coordinate system; p is the 
viscosity; k is the thermal conductivity; and fb represents 
body forces. The above equations are supplemented with 
the state equations 

P = (Y-1)p[e - 2 uiui] (7) 
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outer loop 


T= [e - |uiUi]c v (8) 

for a complete solution, in which y is the ratio of specific 
heats and Cy is the specific heat at constant volume, such a 
formulation being valid for a perfect gas. A solution of the 
nonviscous form of equation (2) may be achieved by first 
obtaining a Taylor series expansion of V in time domain. 
The spatial domain, Q, is next discretized by unstructured 
meshes consisting of 3-D tetrahedron elements. Using 

A A 

linear, finite element approximations, V = aV, V being 
nodal variable values, and employing a Galerkin weighted- 
residual procedure, a time-dependent form of the governing 
equations may be obtained as 

MAV = -At[cM + K]V + R (9) 

in which c is a scalar, and R includes artificial viscosity 
effects essential for capturing shocks. 

Fq^ti nn (9) is solved by advancing the time-dependent 
form until steady conditions are obtained. Both an explicit- 
timestepping, iterative scheme^ and a quasi-implicit 
solution algorithm have been incorporated in the STARS 
program to that effect Thus, in the latter procedure, at a 
typical n+1 1 ** iteration stage, assuming the following 
relations, 

AV = V n+ i - V n 

V = (V n +i - Vn) / 2 
equation (9) may be expressed as 

[(i + yc)M +^k]v„ +1 = 

[(l - y c)M - y K ]v„ + R (10) 

or 

[A]V n+1 = [B]V n + R (11) 

Expressing 

A = D + O (12) 

in which D and O are respective matrices containing 
diagonal and off-diagonal terms of A, the iterative solution 
procedure may be carried out as 


[D]V n+ l = [B]V„ - [0]V n+1 + R (13) 
Solve V n+ i iteratively 
inner loop 

V£V = [D]"*([B]V n - [OlV^ + r) (14) 

If V<£\> * V n +\, continue iteration 
STOP 

If V n+ 1 * v n , go to n+2 step 
STOP 

Both solution schemes prove to be suitable for effective 
solution of practical problems. 

The above solution techniques have been applied to a 
fluids domain idealized by unstructured grids. An 
advancing front technique, developed for automated 3-D 
mesh generation,^ has been found to be suitable for 
discretization of complex domains. This procedure has 
been further modified 3 to effect significant savings in 
solution time; this is achieved by first generating a grid 
whose cells have linear dimension about twice the required 
size, and then reducing each cell locally to reach desired cell 
sizes. 

The STARS multidisciplinary analysis module 
performs aeroelastic and ASE analyses. For linear 
systems, 4 the code uses panel methods to generate unsteady 
aerodynamic forces, and subsequent flutter analyses may be 
performed by either the k, p-k, or ASE method. A far 
more elaborate effort is required, however, for the nonlinear 
case in which the CFD technique yields the unsteady 
forces. Notable earlier efforts in this area have been based 
on the finite difference method, 5 ' 7 whereas our approach 
has been based entirely on the finite element technique. 
Figure 2 shows a flowchart of the major numerical 
solution steps adopted in the STARS program to perform 
nonlinear flutter analyses for complex aerospacecraft. 
References 3 and 8 give in detail related numerical 
formulations in which the generalized equations of motion 
are cast in state-space form from which a response solution 
can be achieved easily by using a standard procedure. A 
fast Fourier transform produces such response solution data 
to obtain the damping value, which is indicative of the 
degree of aeroelastic instability. 

A 

The generalized aerodynamic force vector, f (fig. 2), is 

assembled from nodal pressure values of the finite element 
structural grid which are in turn computed by interpolation 
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from such values at aerodynamic grid points derived from a 
CFD Euler solution. In this process, for each structural 
node pertaining to an element, a triangular aerodynamic 
element encircling the node is first identified and its average 
nodal pressure value assigned to the node. This process is 
repeated for the rest of the nodes of the element and 
averaged among the nodes; such calculations are performed 
for all structural elements. 



Fig. 2 Flowchart for nonlinear flutter analysis. 

In the area of heat transfer analysis, the associated 
module can solve heat conduction problems for a wide array 
of structures, including composite ones. The program is 
capable of handling both linear and nonlinear radiation 
boundary conditions. 


Panel Flutter Analysis 


Reference 9 presents details of extensive flutter 
calculations for a rectangular panel employing approximate 
aerodynamic theory, and such results are also compared 
with experimental findings. The STARS program has 
been used to perform nonlinear flutter analysis of the panel 
with clamped edges, using the CFD and ASE modules of 
the code. Figure 3 shows the surface grid of the 
aerodynamic domain around a plate of aspect ratio 2. 
Figure 4 depicts a comparison of flutter solutions by the 
various procedures; the flutter parameter is defined as X = 


2qa 3 

pD 


, in which 


q = airstream dynamic pressure 


P 

D 


= Vm 2 - 1 , M being the Mach number 
Et^ 

= TT7: — 57* E, t, and v being the elastic 

12(l-v z ) 

modulus, plate thickness, and Poisson's 
ratio, respectively 



Fig. 3 Surface grid for aerodynamic domain for the 
clamped plate problem. 


NUMERICAL EXAMPLES 


Generic Hypersonic Vehicle Aeroelastic Analysis 


The STARS computer program is currently used to 
solve a large number of project-related problems at the 
NASA Dryden Flight Research Facility. Also, a large 
number of test cases has been solved by the program to 
assess the efficacy of solution algorithms and tools. Some 
relevant examples are presented here. 


A new finite element structural model of the vehicle, 
upgraded from an earlier model, 10 was recently generated in 
which special emphasis was placed on the generation of 
well-conditioned elements; the relevant numerical model 
has the following details: 

No. of elements = 4,990 

No. of nodes = 2,812 

No. of degrees of freedom = 16,872 
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Fig. 4 Clamped rectangular panel - comparison of 
experimental, approximate aerodynamic theory and STARS 
nonlinear aeroelastic solution. 

A free vibration analysis of the vehicle was performed by 
the STARS program using a block Lanczos algorithm; 
table I presents a summary of such results, and figure 5 
shows some of the mode shapes. 

Table I. Generic Hypersonic Vehicle - 
STARS free vibration analysis. 


Mote# 

Frequency (Hz> 

Description 


STARS 

LaRC/EAL 


1-6 

0.0 


rigid body modes 

7 

3.18 

3.01 

fus. lb, symm 

8 

4.03 


wing lb, anti 

9 

4.05 

4.02 

wing lb, symm 

10 

5.56 


fus. lb, anti 

13 

6.91 

7.70 

wing It, symm 

14 

6.97 


wing It, anti 

15 

7.27 

7.06 

fus. 2b, symm 

16 

9.40 

9.47 

fin lb, symm 

17 

9.45 


fin lb, anti 


Unsteady aerodynamic forces were next generated, 
using the linear aerodynamic module of STARS, 
employing panel methods. Figure 6 depicts the proposed 
flight envelope of the vehicle, also providing data on flutter 
analysis points, whereas figure 7 shows the vehicle linear 
aerodynamic paneling. Tables II and III give details of such 
analysis results. 

Table II. Generic Hypersonic Vehicle - STARS 
symmetric half-aircraft flutter analysis. 


Altitude (ft) 

Mach No. 

Sneed (TCEAS1 

Freq. fH z) 

Sea level 

0.9 

1807 

5.5 


0.9 

2273 

3.3 

10,000 

0.9 

1798 

5.5 


0.9 

2265 

3.3 

30,000 

2.2 

2286 

5.4 

50,000 

2.2 

2298 

5.3 


3.4 

2693 

5.4 

70,000 

5.6 

3347 

5.6 

Table III. 

Generic Hypersonic Vehicle 

; - STARS 


antisymmetric half-aircraft flutter analysis. 


Altitude (ft) 

Mach No. 

Speed (KEAS1 

Freo. (Hz) 

Sea level 

0.9 

1435 

6.0 


0.9 

1636 

5.5 

10,000 

0.9 

1428 

6.0 


0.9 

1641 

5.5 

30,000 

2.2 

1870 

5.6 

2.2 

2314 

5.2 

50,000 

2.2 

1930 

5.6 


2.2 

2314 

5.2 


3.4 

2271 

5.6 


3.4 

2700 

5.3 


In an attempt to obtain nonlinear flutter solutions in 
the hypersonic flow regime, a set of Euler analyses were 
performed using the STARS CFD module. Thus, figure 8 
shows some views of the associated unstructured 
aerodynamic grid that has the following details: 

No. of elements = 262,787 

No. of nodes = 46,797 

The Euler analyses were performed for a number of flight 
conditions, and figures 9 and 10 show vehicle speed and 
pressure distribution for a typical flight condition, Mach = 
5.6 and a = 0. 
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Fig. 5 Generic Hypersonic Vehicle - some 
typical mode shapes 

(symmetric (S) and antisymmetric (A/S)). 
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(a) Vehicle surface aerodynamic mesh. 



(b) Surface grid for 3-D solution domain. 

Fig. 8 Generic Hypersonic Vehicle - complete surface mesh generation. 
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Fig. 9 Generic Hypersonic Vehicle (Mach = 5.6) - plan view of Mach number distribution. 
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Fig. 10 Generic Hypersonic Vehicle (Mach = 7.0) - plan view of pressure distribution. 
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CONCLUDING REMARKS 

A fmite-element-based analysis procedure developed for 
effective modeling and simulation of aerospacecraft that 
exhibit multidisciplinary interaction has been presented. 
Some numerical examples are also presented that 
demonstrate the applicability of the STARS computer 
program for the analysis of practical problems. Since both 
structural and fluids modeling are accomplished by the 
common finite element method, transfer of data between 
the two systems proves to be a natural process. Although 
a nonlinear aeroelastic analysis for a complete vehicle 
requires rather extensive computing resources, utilization of 
the current finite element-based procedure is justified in 
view of its ability to accurately model complex vehicle 
geometry. 
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